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XII. On the expansion in a series of the attraction of a Spheroid. 
By James Ivory, M. A. F. R. S. 

Read January 17, 1822. 

1 he purpose of this paper is to make some observations on 
the developement of the attractions of spheroids, and on the 
differential equation that takes place at their surface. 

1 . The whole of this doctrine depends on one fundamental 
proposition. Let/ (6, <p) denote any function of the sines and 
cosines of the variable arcs and <p ; and put [a = cos 8 : then 
the given function may be developed in a series, viz. 

/(fl, 9 ) = Q° + Q (,) + Q (2) . • • + Q (i) . • • &c. 
every term of which will separately satisfy this equation in 
partial fluxions, viz. 



•{<—>#} 



«fiH+f(,- + ,)Q» = .. 



dp ' 1 — /** ' dip 1 

Now in one case there is no difficulty. Whenever/(0, <p) 
stands for a rational and integral function of (x, VT^j? . sin q>, 
Vi— . [a 2 , cos <p ; or of three rectangular co-ordinates of a point 
in the surface of a sphere; the proposition is clear. In this 
case the same combinations of the variable quantities are 
found in the terms of the series and in the given function ; 
and by employing the method of indeterminate coefficients, 
the two expressions may be made to coincide. The inquiry 
is therefore reduced to examine the nature of the develope- 
ment when/ ( 6, <p) is not such a function as has been men- 
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tioned. One thing is indisputable. Since the terms of the deve- 
lopement contain no quantities, except such as arise from com- 
bining three rectangular co-ordinates, it follows, that when 
f(&, <p) is not explicitly a function of ^ Vi — yJ* . sin <p, 
Vi — [a\ cos <p, it must be considered as transformed into such 
a function. Algebraically speakings the transformation is no 
doubt always possible ; but there may be danger that, by 
proceeding in this way, we fall upon expressions which are 
not proper representatives of the given function ; which are 
symbolical merely, and which cannot be safely employed in 
the investigation of truth. 

In order to fix the imagination; and to avoid every sort of 
uncertainty and obscurity, I shall take a particular, although 
a very extensive* case of the general expression. IshalL sup- 
pose that/(0, <p), or more shortly y, denotes a rational and 
integral and finite function of the four quantities, sin 6, cos 9, 
sin <z>, cos <p. We may then substitute for the powers and pro- 
ducts of sin q> and cos <p, their values in the.sjnes and cosines 
of the multiples of the arc;- .by which means we shall obtain, 

y = M (o) + M (I) cos <p + M (z) cos 2<p + &c. ; 
+ N^ 1 ' sin <p -f- N (2) sin 2<p 
the symbols M^ o) , M^\ N^ &c. standing for rational ;and 
integral functions of cos Q and sin 6, or of ja and Vi— p/-. Again, 
every even power of Vi — \i? is an integral function of p; 
and every odd power is equal to a similar function multiplied 
by Vi — ^ : the value of y will therefore be thus expressed, 

+ (G ^)+Vi-^.^ (1 ( V))sin <p 
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all the functional quantities being integral expressions of p>. 
The general term of this expression is 

{ F%) + VTZ^./U ] cos i<p 

+ { G C V) + ^1— V>-g % } sin i<p ; 
> 

and by multiplying byii=£Il, or i, it- will become, 

\ 1 ■" E* f ( lT -V ) cos ^ 

I X"" EI K "^ ' sin Zip; 

and finally;, by expanding the denominators, we get, 

M w (i -lA-y^cosV? + N (i) (l-V)"*" sin z<p; 
the symbols M® and N® denoting rational series of the 
powers of [a. By performing these operations in air the 
terms containing <p, and likewise by expanding the radical 
V i—ft m the first term, the value of y will be thus ex- 
pressed ; 

... 1 2' 

> y=M (0) + MP } ( 3—^) 2 cos 0+,M (2) (i-V) 2 cos sp + &c. 

_|_N W (i__ ^) 2 sin <f)+N (2) (l— [x 2 ; 2 sin 39 
all the functional symbols standing for series of the" powers 
of (ju The given expression now consists entirely of combi- 
nations of the quantities [*, V\ — ^ . cos <p, "^l— (x 2 . sin <p; 
that is, it is a function of three rectangular co-ordinates. 

The same end might have been accomplished by a shorter 
and more simple process, which will apply to every function 
of two variable arcs. 
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Put u=V\ — p* . sin <p, %=V\— y? . cos <p ; then 

COS ^ ===== (JLr 

sin 6 = Vi — jT" 

II: 

sm $ = 



V 

COS(p = 



1-^ 

z 



and by substituting these values and expanding the radical 
quantities wherever they occur, y will be transformed into a 
function of \b, u, z, or of three rectangular co-ordinates. The 
former process has been chosen under the idea, that it exhi- 
bits more clearly the quantities which are in a manner extra- 
neous to the function, and are introduced merely for the 
purpose of making it put on a certain form. 

Having now reduced y to a function of three rectangular co- 
ordinates, the developement in question will be obtained by 
the method of indeterminate coefficients already mentioned. 
Nothing more is necessary than to form the quantities 
Q (0) , Q w » Q (2; , &c., giving to each the most general expres- 
sion, and leaving the coefficients indeterminate ; then the 
sum of these terms will contain the same combinations of 
p, Vi — p' . sin <p, Vi — \x. cos <p, thatjy does ; and by making 
the two expressions coincide, we shall obtain,* 

Now let us go back to the original expression of y. By 
following the process described in the second chapter of the 
third book of the Mecanique Celeste, a similar developement of 
that quantity will be obtained. But it is proved, in the same 
chapter, that a given function cannot be developed two diffe- 

* Mecanique Celeste, Tom. ii. p. 42. 
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rent ways, in a series of terms that satisfy the general equa- 
tion in partial fluxions :* and from this it appears, that the 
developement obtained by the algebraic operations described 
above, is identical with the developement found by Laplace's 
method. 

There is another way of expressing the terms of the deve- 
lopement of y, namely, by definite integrals. But upon this 
head there is no difficulty, when the possibility of the deve- 
lopement is allowed. The question is not concerning the pro- 
perties of a series of quantities that satisfy the general equa- 
tion of partial fluxions ; but whether the developement can 
be admitted in all cases as a fit instrument for the inves- 
tigation of truth. 

Since it has been shown that the developement of y, ob- 
tained by the procedure in the second chapter of the third 
book of the Mecanique Celeste, is the same with the like deve- 
lopement found by immediately transforming the given ex- 
pression into a function of three rectangular co-ordinates; it 
is just to say of the former method, that in reality it is nothing 
more than a particular way of effecting such a transforma- 
tion, while at the same time it gives to the transformed quan- 
tity a certain arrangement. In the process we have followed, 
there is no need to employ the differential equation that takes 
place at the surface of the spheroid ; and by thus going more 
directly to the foundations of the method, we can discern 
more clearly what goes on under the cover of many compli- 
cated analytical operations. There is an essential distinction 
with regard to the developement to be observed between the 
two cases when the given expression y is explicitly a function 

* Mec. Cel. Tom. II. pp. 32, and 33. 
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of three rectangular coordinates, and When it must be made 
to assume tM'fofm of such a function by a transformation. 

When y is an explicit function of three rectangular co-or- 
dinates; these things are true : 

ist. The developement contains no quantities except what 
are found injy. The two expressions are entirely equivalent, 
being in reality the same quantities differently arranged. 

sdly. Whenjy is a finite expression, the developement will 
consist Of a finite number of terms ; and when y is a converg- 
ing series of an infinite number of terms, the developement 
will be a like converging series. For, in the case of a con- 
verging series, we may approach to the value of y as near as 
we please, by taking in a determinate number of the terms; 
and the developement of this portion of the series will like- 
wise consist of a finite number of terms. 

On the other hand, the properties that demand attention 
are very different, when we suppose that y is not explicitly a 
function of three rectangular co-ordinates. 

ist. The developement will always contain an infinite num- 
ber of terms. 

sdly. In most instances ; and, more particularly, in the 
very general example that has been considered above; the 
terms of the developement will involve an infinite number of 
quantities which do not appear in the original function, and 
which are introduced merely in order to give the develope- 
ment its peculiar form. 

The original function and the whole infinite series of which 
the developement consists, may be represented by this equa- 
tion, viz. 

y = y -f- M — M : 
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the right hand side, which stands in place of the developement, 
containing all the quantities in y, and besides an infinite num- 
ber of other quantities, denoted by -f- M and — M, which 
have opposite signs, and destroy one another. But this mutual 
destruction of the extraneous quantities will take place only 
when the totality of the series, comprehending its infinite 
number of terms, is taken into account. Any determinate 
number of terms of the developement will contain quantities 
not to be found in y ; and, for aught that appears, the diffe- 
rence may have any given amount. No finite part of the 
developement can therefore properly be said to represent the 
proposed function. If we take a separate term, as Q (i) , it is 
not extravagant to say that it may have nothing in common 
with the original expression. Nevertheless such a term is a 
necessary part of the series, in order to balance quantities 
that occur with opposite signs in other terms. 

For the sake of illustration, suppose a spheroid of revolu- 
tion determined by this equation, viz. 

e being a small coefficient, and ^ the cosine of an arc reckoned 
from one of the poles of revolution. In such a spheroid the 
greatest diameters will make an angle of 45 with the equa- 
tor, resembling the figure which the planet Saturn was, some 
years ago, supposed to have. In order to develope the vari- 
able part of the radius of the spheroid, that is |*Vi— p,', in a 
series of quantities that satisfy the equation in partial fluxions, 
it is first necessary to expand Vi—^; and the number of 
terms of the developement will therefore be infinite. 

To the preceding equation let there be added the term 

mdcccxxii. p 
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fy. sin <p ; <p denoting the variable angle that the circle on 
which the distance from the pole is reckoned, makes with a 
circle given by position; and the nature of the spheroid will 
be thus expressed, viz. 

rsszall -f- e(yVi — p +//* sin <p) X. 

In this case the quantity to be developed must be put under 
this form, viz. 

and the developement will not only consist of an infinite num- 
ber of terms, but these terms will contain an infinite number 
of quantities which arise from the expansion of the radical in 
the denominator, and which are not to be found in the ori- 
ginal function. 

There is therefore a real distinction to be made between 
the two cases when y is an explicit function of three rectan- 
gular co-ordinates, and when it is not. A method of calcula- 
tion which is clear, exact and elegant, when it is confined to 
the first case, becomes clouded with obscurity, if not merely 
symbolical, when it is extended to the other case. To say 
the least, there are certainly great difficulties which are not 
explained ; and if there be any geometers who hesitate, and 
have doubts, they are not without their excuse, and ought 
not to be entirely condemned. 

2. We come next to consider the differential equation that 
takes place at the surface of a spheroid. Of this equation, 
three demonstrations have been published ; one, in the second 
chapter of the third book of the Mecanique Celeste ;* another 
by the same author, not precisely the same with the former, 

* Tom. II. p. 28. 
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but similar to it, in a memoir read to the Academy of Sci- 
ences in 1818 ; and a third by M. Poisson, in an interesting 
and profound memoir on the distribution Of heat in solid 
bodies. The two last demonstrations are fundamentally the 
same ; but as M. Poisson has stated the reasoning more 
fully, and fixed the sense of the proof more precisely, I wish 
to refer to his memoir. One observation it is proper to make, 
which is, that in the integration by which the differential 
equation is proved, the function expressing the thickness of 
the molecule is considered as invariable, or is treated as a 
constant quantity. It is essential to attend to this remark, 
which in reality affords the clue necessary to unravel what 
is mysterious in this investigation. 

In order to acquire a distinct notion of the meaning of the 
differential equation in the sense in which it is demonstrated, 
conceive the surface of the earth, perfectly smooth and spheri- 
cal, to be covered with circles, we shall say, of a thousand 
yards radius each. The circles may either touch one ano- 
ther and cover the whole surface of the earth ; or they may 
cover it partially only, and with any interruptions that can 
be imagined : conceive also that a mass of matter, or mole- 
cule, is placed within every circle ; the thicknesses of the 
molecules being entirely arbitrary, and subject to no law of 
variation or restriction whatever, excepting that they are 
quantities of inconsiderable magnitude when compared with 
the radius of the sphere. These things being supposed, the 
differential equation will be separately true of every one of 
the molecules. 

Let now the whole surface of the earth, or any portion of 
it, be covered with molecules, the thickness varying according 
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to a particular law, or being expressed by some function of the 
arcs that determine the position of a molecule : the differen- 
tial equation will be true of all these molecules. But it will 
be equally true of them, if the thickness be entirely arbitrary, 
and subject to no law of variation. It is true, in reality, of 
each molecule taken in an insulated manner, and because its 
thickness is some determinate quantity. How then can such 
an equation be a fit means of proving that the thickness 
varies in one way rather than another, or that it comes under 
a particular developement ? 

To consider this matter more particularly, let y =/(0, <f>) 
denote, as before, the thickness of the molecule ; and put, 

y = cos cos 6' -J- sin 6 sin 6' cos (<p — <?>'), 

f = Vtf — 2ra . y + a% 
a being less than r: then tt denoting the semi-circumference 
to the radius l, the differential equation at the surface is what 
the following formula becomes in the particular case of r=a, 
viz. 

d 'l\ 
or, by substituting the value of ^ , 

\ r* — 2ray + a* } 2 

the integral being taken between the limits 6'=o, <p'=o and 
0':=7r, ^'==2fl-, and making a=r, after the integration. 

Now, when r=a, r* — a 2 =o; and in this case all the ele- 
ments of the integral are evanescent, unless in the particular 
circumstances when the denominator is evanescent, or when 



in a series of the attraction of a spheroid. 109 

iy=i ; which can happen only when 6' =0, and <?>'=$>. If 
therefore we suppose r* — a 2 to be indefinitely small, the 
whole value of the integral will be obtained by extending the 
integration to such values of 0' and <p' as differ indefinitely 
little from 6 and <p. But when & and <p' differ indefinitely little 
from and <p,y' orf(Q', <p'), will differ indefinitely little from 
y, or/ (5, <p) : and hence it follows that, in the whole extent 
of the integral, we may consider y' as constant and equal to y. 
We have therefore to prove the truth of this formula, viz. 

_ yr /Y*(»* — fll ) s in 8' d& dp 
y — ^'JJ ~, \1 ' 

\r»— 2ray+o*/ z 

between the limits 0'=o, <?>'=o, and (?=#, <j>'=27r, in the par- 
ticular case of a—r. 

The two arcs 0, 0' and the arc of which y is the cosine, are 
the three sides of a spherical triangle ; <p— 9' is the angle op- 
posite to the last arc ; and if ^— \|/ denote the angle oppo- 
site to 0', the element of the spherical surface will be equally 
expressed by sin Q'dQ'dtp', or dtfd. cos 0, and by di/dy. where- 
fore we have to integrate this formula, viz. 

yr pr* (r 1 — fl a ) . d\'dy 

y—^jj 1 u 

\y*— aroy+a*/ 2 

Now integrate between the limits \{/=o and \|/= 27r : then 

v — E. C fr*— a») dy m 

y~2.ji \a • . 

[ri—zray+a*) 2 

integrate again, and 

y r* — a 7 - 

* za ' \/Z 



Vr* — zray+a % 

then take this integral between the limits y = + 1, and 
<y = — 1 ; and 

y_ f r % —a % r*—a % \ 

' za ' \ r — a r+« J' 
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which equation is true when we make r = a. This analysis 
is equivalent to the demonstration of M. Poisson. Whatever 
may be thought of the reasoning, it cannot be denied that, in 
both processes, y' is in fact treated as a constant quantity. 
The equation is true of each individual molecule taken sepa- 
rately, and merely because its thickness has some determinate 
value. Such a demonstration cannot therefore be employed 
to prove that the thickness of a series of molecules covering 
the surface of a sphere, or a part of that surface, follows a 
certain law of variation, or comes under a particular deve- 
lopement. r... -<-.-.,. 

But a legitimate process of reasoning requires that, in the 
formula (A), while a represents any determinate quantity 
less than r, y' be considered as a function of the variable 
quantities sin 0', cos 6', sin <p', cos <p' ; and likewise that the 
integration be extended to the whole surface of the sphere, 
or to that part of it covered with the related molecules ; after 
which the true value of the formula will be obtained by 
making a — r. The whole system of molecules being com- 
prehended in the result; we may thence deduce, by a reverse 
process, the law according to which their thickness must vary, 
in order to produce that result. Now the integration here 
spoken of, cannot be executed, unless in the case when y' is 
explicitly a function of three rectangular co-ordinates. It is 
therefore only in this case that the differential equation can be 
considered as rigorously , proved ; and it is remarkable, that, 
when we seek from that equation the developement of y', it 
always comes out in a function of three rectangular co- 
ordinates. 

Wheny is not explicitly a function of three rectangular 
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co-ordinates, the formula (A) cannot be integrated. And, 
perhaps, what is now said, is alone sufficient to show that, in 
this case, some modification takes place, which it were de- 
sirable to have fully explained. On attempting to transform 
y into an expression containing <y, \/i — y % , sin tf/, cos -A/ in 
place of cos 0', sin 0', sin <p', cos <p', the powers of 1/1— y* 
make their appearance as divisors ; and hence it is to be feared 
that the integral will be infinite at the limits ; which circum- 
stance would make it impossible to conclude with certainty 
what the value sought will become in the particular case of 
a = r. But it would be of no utility to seek a strict demon- 
stration of the differential equation : because in reality the 
method, when it is extended to all functions of two variable 
arcs, is independent of that equation, being derived from this 
proposition, that every such expression is either explicitly a 
function of three rectangular co-ordinates, or may be trans- 
formed into one. The developement in question may always 
be found, as has been shown, by the rules of algebra ; and 
the differential equation is wanted neither for proving the 
possibility of the developement, nor for calculating its terms. 
But in this plainer way of considering the matter, it appears 
that the developement does not represent the given expression 
y, when that expression is not an explicit function of three 
rectangular co-ordinates, in the same sense that it does when 
it is such a function. There is, therefore, a difficulty left 
unexplained ; and we may be permitted to doubt, whether so 
important a part of the celestial mechanics, as that regarding 
the figure of the planets, rests, with sufficient evidence, on 
the doctrine laid down concerning the generality of the 
developement. 
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The observations that have been made, relate only to the 
foundation of the method in the Mecanique Celeste, and do 
not touch upon the general scope of the analysis, which is 
deserving of every praise. It is natural to think, that the 
theory of the figure of the planets would be placed on a 
firmer basis, if it were deduced directly from the general 
principles of the case, than when it is made to depend on a 
nice, and somewhat uncertain point of analysis. To hazard 
a conjecture suggested in the course of writing this paper, 
the theory will probably be found to hinge on this propo- 
sition, that a spheroid, whether homogeneous or heterogene- 
ous, cannot be in equilibrium by means of a rotatory motion 
about an axis, and the joint effect of the attraction of its own 
particles, and of the other bodies of the system, unless its 
radius be a function of three rectangular co-ordinates. If 
this proposition were clearly and rigorously demonstrated, 
the analysis of Laplace, in changing the ground on which it 
is built, would require little or no alteration in other re- 
spects. 

Dec. 22, :82i. 



